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Superfluidity and superconductivity are genuine many-body manifestations of quantum coherence. 
For finite-size systems the associated pairing gap fluctuates as a function of size or shape. We provide 
a theoretical description of the zero temperature pairing fluctuations in the weak-coupling BCS limit 
of mesoscopic systems characterized by order/chaos dynamics. The theory accurately describes 
experimental observations of nuclear superfluidity (regular system), predicts universal fluctuations 
of superconductivity in small chaotic metallic grains, and provides a global analysis in ultracold 
Fermi gases. 
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A microscopic theory of superconductivity based on 
pairing was set up in 1957 by Bardeen, Cooper and Schri- 
effer These theoretical ideas were subsequently ap- 
plied to finite systems by Bohr, Mottelson and Pines to 
describe ground-state superfluid properties of atomic nu- 
clei Q . Today pairing effects are central in a broad range 
of quantum systems, including neutron stars, metallic 
grains, atomic gases, nuclei, etc 0, 0, H, H|. As the sys- 
tem size diminishes, finite-size effects become important 
and lead to corrections with respect to the bulk homo- 
geneous behavior. Of particular interest is the influence 
of the discreteness of the single-particle quantum energy 
levels. In connection with superconductivity, its impor- 
tance was initially emphasized by P. W. Anderson Q, 
who pointed out that superconductivity in small metal- 
lic grains should disappear when the single-particle mean 
level spacing becomes of the order of the pairing gap. The 
validity of this criterion was qualitatively confirmed ex- 
perimentally in the 90's 0]. Another consequence of the 
discreteness of the energy levels is the appearance of fluc- 
tuations as a parameter is varied. There is at least one 
clear experimental evidence of fluctuations of the pairing 
gap in superfluid systems, through the odd-even stagger- 
ing of nuclear masses as a function of the mass number. 

Our purpose here is to present, within a mean-field 
approximation, a theory of the pairing-gap fluctuations 
valid in the weak coupling BCS limit for arbitrary ballis- 
tic potentials. Our method is based on periodic orbit the- 
ory, which has been successful in describing mesoscopic 
fluctuations of thermodynamic and transport properties 
in many-body systems Q. The results allow for a de- 
tailed calculation of the fluctuations in specific systems. 
In particular, they provide an accurate description of the 
pairing fluctuations in nuclei (cf Fig. [5]). We also focus 
on statistical properties, which are shown, generically, 
to be non-universal. The analysis leads to a global and 
complete picture of the typical size of the fluctuations in 
terms of properties of the corresponding classical system, 
namely regular or chaotic dynamics. 

Our starting point is the mean field BCS equation for 



the pairing gap A [lj], 
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where G fixes the strength of the pairing (seniority) in- 
teraction, p(e) is the single-particle level density, and we 
have put the Fermi energy to zero. The energy cut off ±L 
is given by the physical conditions, that are often related 
to the determination of the force strength G. Following 
semiclassical approaches, we divide the pairing gap as 
well as the single-particle density of states in a smooth 
part and a fluctuating part, A = A + A and p = p + p, 
respectively. In the weak coupling limit A <C L, the 
smooth part of the gap is given by the well known solu- 
tion A = 2Lexp(— 1/pG) (see Ref. [lOj for regularization 
schemes). The fluctuating part of the density p can be 
expressed as @ 

oo 

p(e) = A p , r cos(rS p /h + v Pi r), (2) 
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where the sum is over all primitive periodic orbits p 
(and their repetitions r) of the classical underlying ef- 
fective single-particle Hamiltonian. Each orbit is charac- 
terized by its action S p , stability amplitude A p ^ r , period 
t p = dSp/de and Maslov index v p r (all evaluated at 
energy e). Assuming A <C A, an equation for A may 
be obtained by multiplying Eq. ([T]) by A, replacing the 
single-particle level density by its semiclassical expres- 
sion, and expanding up to lowest order in fluctuating 
properties. Assuming moreover Ac L gives 
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This equation, where all classical quantities involved are 
evaluated at Fermi energy, contains detailed information 
about the variations of the pairing gap. Note that A only 
depends on G and L through A. K (x) is the modified 
Bessel function of second kind. Through it, a new char- 
acteristic "pairing time" associated with the pairing gap 
is introduced, 



(5) 
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Since Kq(x) oc exp(— x)/y/x for x 3> 1, the Bessel func- 
tion exponentially suppresses all contributions for times 
t ^> t a (making the sum convergent). The average 
part of the gap A is thus playing, in this respect, a 
role very similar to the temperature in the general the- 
ory of mesoscopic fluctuations (cf Ref. (HI). ^ n contrast, 
K (x) oc — log(x) for and short orbits (compared 

to t a ) are logarithmically enhanced. 

Since the value of the actions depend on the shape 
of the mean-field potential, Eq. ([3]) predicts, generically, 
fluctuations of the pairing gap as one varies, for instance, 
the particle number, or the shape of the system at fixed 
particle number. The fluctuations result from the inter- 
ference between the different oscillatory terms that con- 
tribute to A. The symmetries of the potential and the 
nature (integrable or chaotic) of the underlying classical 
motion are crucial to understand the interference pat- 
tern. When the motion is regular (integrable), contin- 
uous families of periodic orbits having the same action, 
amplitude, etc, exist. The coherent contribution to the 
sum (J3j) of these families of periodic orbits produces large 
fluctuations. In contrast, in the absence of regularity or 
symmetries, incoherent contributions of smaller ampli- 
tude coming from isolated unstable orbits are expected. 
Moreover, aside the dependence on the regular or chaotic 
nature of the single-particle motion, the presence or ab- 
sence of universality in the statistical properties of the 
fluctuations will depend on the dominance of short or 
long periodic orbits. 

We will make below an analysis of the predictions of 
Eq. Q in the nuclear the neutron number is 

varied. Before, and in order to avoid at this stage a de- 
tailed study of a particular system, we concentrate on 
a global analysis, namely the typical size or root mean 
square (RMS) of the BCS gap fluctuations in a generic 
mesoscopic system. The second moment of the fluctua- 
tions may be expressed from Eq. ([3]) as 
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where r H = h/Sis Heisenberg time (6 — p _1 is the single- 
particle mean level spacing at Fermi energy), and K{t) 
is the spectral form factor, i.e. the Fourier transform of 
the two-point density-density correlation function [13]. 

The structure of the form factor K(t) is characterized 
by two different time scales. The first one, the smallest 



of the system, is the period r m i n of the shortest periodic 
orbit. The form factor is zero for r < T m ; n , and dis- 
plays non-universal (system dependent) features at times 
T min ~ t < t h . As t further increases, the function be- 
comes universal, depending only on the regular or chaotic 
nature of the dynamics, and finally tends to t h when 
r 3> t h . The result of the integral (O thus depends on 
the nature of the dynamics, and on the relative value of 
t a with respect to r m j n and t h . According to Ander- 
son criterion superconductivity exists if A > S (we 
are not interested here in the ultrasmall regime A < 6 
[r| where the BCS theory fails). Then, A > S implies 
Ta/t h = 6/2irA -C 1. Because the Bessel function Kq ex- 
ponentially suppresses the amplitude for times r 3> t a , 
one can safely ignore the structure of the form factor for 
times of the order or bigger than Th, and use the so called 
diagonal approximation of K(t) 12]. In the simplest ap- 
proximation, all the non-universal system-specific fea- 
tures are taken into account only through r m i n [Til ] . and 
one can write K(t) = for r < r m i n and, for r > r m i n , 
K(t) = t h for integrable systems and K(r) — 2r for 
chaotic ones with time reversal symmetry. 

This finally gives the expressions for fluctuations of 
the pairing gap (normalized to the single-particle mean 
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level spacing), a = < (A 2 ) /5, assuming regular dynam- 
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where .F„(L>) = 1- J° x n Kg(x)dx/ J™ x n Kg(x)dx. The 
argument 
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is a system dependent quantity inversely proportional 
to the dimensionless conductance, g = T H /r min , an in- 
trinsic characteristic of the system independent of the 
pairing coupling. D can also be viewed as the system 
size divided by the coherence length of the Cooper pair, 
£o = JiVp/(2A), where v F is Fermi velocity. Equations 
(|7|) and (|H1) , which together with Eq. J3]) are the main re- 
sults of this study, show that the variance of the pairing 
gap is a function of its normalized mean part, A/<5, and 
of the dimensionless conductance, g, as shown in Fig.[TJ 
The monotonic function F n (D) has the following limit- 
ing behaviors, F n {D) — > 1 as D — > 0, whereas F n (D) — > 
exponentially fast for D 3> 1. Thus, in a system charac- 
terized by large g- value, D — > 0. In this case F n (D) — ■+ 1, 
all dynamical system specific properties disappear, and 
we obtain a "universal" (system independent) behavior 
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FIG. 1: Fluctuations of the pairing gap as a function of the 
mean value for mesoscopic systems (log-log scale; all quanti- 
ties normalized with the single-particle mean spacing). The 
dashed (regular, in blue) and dotted (chaotic, in red) curves 
correspond to different values of the dimensionless conduc- 
tance, g, and the limiting case of g — > oo is shown by solid 
lines. The results are valid in the pairing regime A/5 > 1. 
Applications to Nuclei and Nano-grains are marked out. Ul- 
tracold atomic gases may be controlled to appear in major 
parts of the figure. 
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FIG. 2: Nuclear pairing gaps for neutrons. The blue dashed 
line shows average experimental gaps, and isotopic variation 
within one standard deviation is marked by the shaded area. 
Calculations from the cavity model are shown by the solid 
red line. Average behavior (Eq. {TO}) is shown by solid black 
fine. Inset: RMS of pairing gap fluctuations versus mass num- 
ber A. The dots are experimental data for protons and neu- 
trons. Solid and dashed lines are regular (Eq. (0) and chaotic 
(Eq. ((8}) predictions, respectively. Data from Ref. [24j ] . 



of the gap fluctuations given by the prefactors in Eqs. ([7|) 
and that correspond to a pure uncorrelated Poisson 
sequence and to a GOE random matrix spectrum, respec- 
tively (the latter, a 2 ch = 1/(2tt 2 ), was obtained previously 
in Ref. [Hj]). This situation is shown by the solid lines 
in Fig. [TJ purely GOE fluctuations imply a constant am- 
plitude of the normalized fluctuations of the pairing gap, 
whereas an increase with A/S is seen for systems with 
uncorrelated spectra. In contrast, in the generic case 
of systems characterized by finite values of g, F n , and 
therefore the pairing fluctuations, may significantly devi- 
ate from universality (cf Fig. [I]). Thus, in general, pure 
statistical models (like GOE) do not provide an adequate 
description of the pairing fluctuations. 

We shall now apply these results to different physical 
situations. Our first example is a system dominated by 
regular dynamics, namely ground states of atomic nu- 
clei, which bring the best experimental data available 
at present on the superfluidity of finite Fermi systems. 
The ground-state superfluidity of atomic nuclei implies a 
mass difference between systems with an even and odd 
number of particles. The connection between the pairing 
gap and the mass differences is given by the three-point 
measure A 3 (M) = B(M) - [B(M + 1) + B(M - l)]/2, 
where M is the odd neutron N or proton Z number. In 
the presence of other possible interactions, this quantity 
has been shown to be a very good measure of pairing cor- 
relations [l5(. A3 is shown in Fig. [2] for neutrons. The 
average dependence of the neutron and proton gaps is 
well approximated, from experimental data, by 
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2.7 



MeV 
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where A = N + Z is the total number of nucleons [16| • 
We notice a rather strong variation around the average 
value. The A dependence of the RMS of the experimental 
pairing fluctuations is shown in the inset of Fig. [5] 

In order to evaluate the RMS of the pairing fluctua- 
tions from the theoretical expressions, Eqs. © and ©, 
we need the following estimates of nuclear properties (for 
one nucleon type): mean level spacing, S ~ 50 A4 MeV, 
and dimensionless conductance g ~ I.6A 2 / 3 [17]. This 
gives D w 0.21A 1 / 12 , which ranges from 0.27 for A = 25 
to 0.33 for A — 250. Though these values clearly set 
atomic nuclei in the regime T m j n /r A < 1 where pair- 
ing fluctuations are important, we are still far from 
D = T m j n /T A = 0, so that significant deviations from 
universality are expected. By inserting these estimates 
in Eqs. ([7]) and ([5]) the fluctuations are easily evaluated, 
assuming regular or chaotic dynamics. The resulting 
curves are compared to the experimental one in the in- 
set of Fig. O Note, as expected [18|, the good agree- 
ment between the regular dynamics and the experimen- 
tal curve, either in the overall amplitude as well as in the 
A-dependence. 

One may go beyond a statistical description, and use 
Eq. ([3]) to obtain a detailed description of the fluctu- 
ations. For that purpose, we assume for the nuclear 
mean field a simple hard-wall cavity potential. The shape 
of the cavity at a given number of nucleons is fixed by 
minimization of the energy against quadrupole, octupole 
and hexadecapole deformations. To simplify, we take 
N = Z. The periodic orbits of a spherical cavity are 
used in Eq. ©, with modulations factors that take into 
account deformations and inelastic scattering [19j| . We 
set the average of A to zero, as was done with the exper- 
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imental data (although of interest by itself, we will not 
consider its behavior here). In Fig. [2] we compare the the- 
oretical result A(N) to the experimental value averaged 
over the different isotopes at a given N. The agreement 
is excellent; the theory describes all the main features 
observed in the experimental curve. 

Our second example are the superconducting proper- 
ties of nano-sized clean (ballistic) metallic grains [20(, 
where we may expect the dynamics to be chaotic, see 
Ref. The existence of a superconducting gap was 
demonstrated in the regime A > 5 Q, whereas no gap 
was observed when A < 6 (the transition occurs around 
N ~ 5000, where N is the number of conduction electrons 
in the grain) . The N dependence of the average gap A is 
poorly understood. We will adopt for grains the thin-film 
value A w 0.38 x 10 ~ 3 eV 0|. The mean level spacing 
is S = (2E P )/(3N) w 2.1/ N eV, whereas g w 2.6N 2 / 3 . 
Eq. flSJ) gives D w 4.4 x 10~ 4 N 1 / 3 , which ranges from 
0.05 to 0.02 when N varies between 10 3 and 10 5 . This 
means that the variance will be close to its "universal" 
value obtained by setting F\(D) = 1 in Eq. (JSj> , namely 
cr^ h = 1/27T 2 . Their typical range of variation is repre- 
sented in Fig. [TJ 

In the case of chaotic dynamics, we don't have explicit 
experimental or numerical data to compare with. We 
can, alternatively, compute the fluctuations of the con- 
densation energy, defined as the total energy difference 
between the paired and unpaired system. In the universal 
chaotic limit, our results are in good agreement with the 
numerical calculations of Sierra et al. [2l| , where they use 
Richardson's solution of the pairing problem and random 
matrix theory (GOE) for generating the single-particle 
spectrum. 

Superfluidity in ultracold atomic gases is currently in- 
tensively studied, and provides our third example. The 
confinement potential can be externally controlled to cre- 
ate regular as well as chaotic dynamics, and the atom- 
atom interaction, a, can be tuned around the Fesh- 
bach resonance. Since both particle number and interac- 
tion strength are experimentally controlled parameters, 
the fluctuations may appear in major parts of Fig. [1] 
We estimate S = (2E P )/(3N) and g = ±(3A0 2/3 ; in 
the dilute BCS region A = {2/e) 7 l 3 E F exp {-n/2k F \a\) 
[22l |. with k F the Fermi wavevector, giving D = 
2 7 r(2/e) 7 ' /3 (37V) 1 / 3 exp(-7r/2fc f \a\). Recent experiments 
using Li 6 reach k F \a\ = 0.8 [23|, implying negligible fluc- 
tuations for typical values of N ~ 10 6 . Reducing to 
fe^al = 0.2 and TV = 10 4 yields for generic regular sys- 
tems fluctuations that are on the same magnitude as the 
mean pairing gap, cr reg « 0.5A/5 (lH ]. 

To conclude, we have presented an explicit semiclas- 
sical theory for the pairing gap fluctuations. These are 
generically dominated by system specific features not in- 
cluded in purely statistical models. Different possible 
regimes, as well as the influence of order/chaos dynam- 



ics, were investigated, in particular for the typical size of 
the fluctuations (Fig. [J). The present theory provides, 
for the first time, analytic predictions, valid for a wide 
range of physical situations; it also compares very favor- 
ably with available experimental data. 
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